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Abstract 

r] A joint characterization of reachability (controllability) and observ- 

ability (constructibility) for linear SISO nonuniformly sampled discrete 
systems is presented. The work generalizes to the nonuniform sampling 
Cu the criterion known for the uniform sampling. Emphasis is on the nonuni- 

P form sampling sequence, which is believed to be an additional element for 

I I analysis and handling of discrete systems. 

^ 1 Introduction 

vr^^ The concepts of controllability and observability, first introduced by Kalman [7] , 

^T^ are still two fundamental questions in modern control theory. On the other hand, 

I * the enormous increase in the use of digital computers has stimulated studies 

^-^ in the field of discrete systems. In particular, the problem of controllability 

^^ (reachability) and observability of discrete systems has been already treated in 

^-H the literature in a generalized form [I], [2], - [10] . 

^ Most of the previous references are only concerned with uniformly sampled 

discrete systems. However, the general case of nonuniform sampling offers a 
wider range of situations in the analysis of these concepts for discrete systems, 
i^ The present note tackles this problem and, in fact, a joint characterization of 

reachability and observability for linear single-input single-output nonuniformly 
sampled discrete systems has been developed. Under several conditions, a right 
choice of the sampling instants would guarantee the aforementioned internal 
properties. 

This note emphasizes the importance of the nonuniform sampling sequence 
against other system parameters. Besides the above considerations, nonuniform 
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sampling is believed to be an alternative solution to relevant problems such 
as propagation of measuring and/or rounding errors, communication delays in 
complex computer-controlled systems, identification, etc. 

2 A Joint Criterion for Reachability and Ob- 
servability of Nonuniformly Sampled Discrete 
Systems 

2.1 Model Description 

Consider a linear time-invariant SISO dynamic system 

X{t) = AX{t) + hu{t) Xo = X{0) 

Y{t) = cX(t) (1) 

where X e M" denotes the state vector and u, y G M are the scalar input and 
output, respectively. A{n x n), b{n x 1), c(l x n) are real and constant matrices 
and n is the order of the system. 

As additional assumption, the realization {A, b, c) is required to be minimal. 

2.2 The Reachability Problem 

Let (A, b, c) be an arbitrary minimal realization of order n for the kind of system 
under study. The solution of the state-space equation at the sampling instant 
i„ can be written as 



X{tn) = exp{Atn)Xo + [Gn-1, . • . , Go] 
where 



Un-l 



Uq 



(2) 



G^ ^ exp{A{tn - h))b (i = 0,...,n-l) (3) 

and Uj is the value of the impulse input at time U. 

In mathematical terms, the condition of n-reachability will be accomplished 
if and only if the matrix [Gn~i, ■ ■ ■ , Gq] has full rank. 

Vectors Gi can be rewritten as 

G^^ Bexp{J{tn-U))yo (4) 

where J is the Jordan canonical form of the matrix A, B is the invertible matrix 
of the change of basis, and 

yo ^B-'b. (5) 

Therefore, in order to guarantee the n-reachability we must compute the 
value of 



with 



det[exp{Jajn)ya] 



C^ni ^n — 1 ^n- 



(to = 0, . . . ,n ~ 1) 

i_i (ao = 0) 



and ensure that such a determinant is nonnull. 
First, we denote the components of j/o by 



(6) 
(7) 

(8) 



where rrij (j — 1, . . . , r) is the multiphcity of the r different eigenvalues of the 
matrix A with r < n. (' denotes the transpose). 

Then, we denote the reachabihty canonical form |S] of this arbitrary minimal 
realization by (Are, &re, c^e)- Finally, the impulse response [B] for this kind of 
system can be written as 



h{t) = cexp{At)b — y ^ Ciipi{t) 



(9) 



i=l 



where q € C are constant coefficients and (/?i;M — >■ C (i = 1, . . . ,n) is the 
fundamental system of solutions of an nth-order homogeneous linear differential 
equation. 

Now, making use of the Laplace's expansion by minors, the determinant of 
([6| can be factorized as follows: 



det[exp{Jam)yo] = NiN2det[(p^{am)] 
Let us consider separately each one of these factors 
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(10) 



(11) 



The term A^i is related to the multiplicity of the eigenvalues of the matrix 
A and will always be nonnull 



N2 ^ det 



y\ 



2/ml 



Vml 



2/1 



(12) 



By similarity transformations 

yo = B-H = B-Xe = iCl,C2, . . . , CnY 



(13) 



where Bre is the matrix of the change the basis of Are to the Jordan canon- 
ical form. The term N2 is related to the weighting coefficients (Ci) of the 
characteristic modes in ^ . Remark that A^2 will be nonnuU if and only if 

VL.^O (j = l,...,r). (14) 

This holds, according to the previous meaning of the components of yoj 
because only minimal realizations are considered. 

Finally, [(^^(am)] (i = 1, . . . , n; to = 0, . . . , n — 1) is a n x n matrix involving 
characteristic modes and sampling instants. From ( |10[ ), the following result is 
derived. 

Lemma 1: An arbitrary minimal realization (A, 6, c) is completely w-reachable 
(reachable in n steps) if and only if n consecutive sampling instants are chosen 
in such a way that 

dei[(/?i(a„)] 7^ (i = 1, . . . ,n; to = 0, . . . ,n - 1). (15) 



Proof: The proof is evident from the factorization of ( 10 ) 



Note that if the scalar input had been defined as a control of the form 

u{t) = u{ti) ^Ui ti <t < ti+i (16) 

which is more realistic, then the above result would still be valid. In fact, 
the presence of a data hold does not affect the system characteristic modes. 

2.3 General Criterion 

By duality, the n-observability characterization is straightforwardly derived. 
Indeed, a similar lemma (concerning the n-observability) can be stated by just 
changing the term reachability to observability in Lemma 1. 

The results obtained previously can be unified as follows. 

Theorem: An arbitrary minimal realization [A, b, c) for the kind of system 
under study is jointly n-reachable and n-observable if and only if n consecutive 
sampling instants, not necessarily equidistant, are chosen in such a way that 

det[ipi{a,n)]^0 (i == 1,.. .,n; TO = 0, .. .,n- 1). (17) 

Proof: The proof is evident from the two previous lemmas. It must be noticed 
that the condition (17) depends exclusively on the system characteristic modes 
and the sampling instants. 

Condition (17) imposes a rather weak restriction for the choice of the sam- 
pling instants. In fact, time intervals can be specified so that complete reach- 
ability and observability are preserved. Particular to the uniform case, (17) 
becomes the condition imposed on the sampling interval T (see [5], [5]) when a 
uniform sampling is considered. 

We can also remark that, for this kind of system, the n-reachability and 
n-observability are inseparable concepts. These systems are either n-reachable 
and 7i-observable or, if not, they are neither n-reachable nor n-observable. 



On the other hand, the «-controllabihty (n-constructibiUty) of an arbi- 
trary minimal reahzation can be considered as a weakening of the condition 
of n-reachabihty (n-observabihty) [6]. Therefore, the pair n-controUability/n- 
constructibihty can be characterized as a corollary of the preceding theorem. 
Indeed, both properties will be accomplished if and only if 

exp(yli„)XoeM[G„_i,...,Go] (18) 

where M[. . .] denotes the range space of the columns Gi. 

From the previous theorem and following a similar reasoning, the next corol- 
lary can be proved. 

Corollary: An arbitrary minimal realization (^4, 6, c) for the kind of system 
under study is jointly n-controllable and n-constructible if and only if n -I- 1 
sampling instants, not necessarily equidistant, are chosen in such a way that 

((/7i(a„), . . . , V5„(a„))' e M[((y9i(a„), . . . , <^„(a„))'] (19) 

(with a„j, ipi defined as before) {m = 0, . . . , n — 1) 

a.n = tn - to. (20) 

Remark that the condition of n-controllability / n-constructibility involves 
one sampling instant more than in the preceding characterization. Note also 
that the pair n-reachability / n-observability implies the pair n-controllability/ 
n-constructibility, but the converse may not be true. 

3 A Simple Strategy for Choosing Sampling In- 
stants 

In order to give more insight on the restriction imposed by (17), a simple ex- 
ample is presented. Consider a 2nd-order realization where the matrix A has a 
pair of complex eigenvalues 

a + jb (6>0). (21) 

Now, (6) can be written as 

det[Ya,Yi] = det[exp{Jao)yo, exp{Jai)yo] (22) 

where 

J=[? "M (23) 



and ai,yo are defined as before. 

A geometric interpretation is very easy. Indeed, the generic operator exp(Ja) 
applied on the vector j/o can be viewed as follows: 

1. a counter-clockwise rotation through ba rad; 
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Figure 1: Relationships among the vectors Yo,Yi,Y2 



(c) 



2. a stretching (or shrinking) of the length of yo by a factor exp{aa) 

Thus, a necessary and suficient condition to preserve complete reachability 
and observability is that 

b{ti-to)^kIl A: = 0,1,... . (24) 

Remark that given a first and arbitrary sampling instant tg, just point values 
on the time axis will be forbidden. Thus, (17) is actually not very restrictive. 

On the other hand, given three successive sampling instants (^0:^1:^2) and 
their corresponding vectors (Yq, Yi, 1^2), the Fig. IT] can be interpreted as follows. 

Case a: The sampling instants are chosen in such a way that the pair reach- 
ability /observability is preserved and also the pair controllability /con- 
structibility . 

Case b: The sampling instants are chosen in such a way that the pair reach- 
ability /observability is not preserved but the pair controllability/ con- 
structibility is preserved. 

Case c: The pair reachability /observability is neither preserved nor the pair 
controllability /constructibility. Remark that this situation will never oc- 
cur for uniform sampling. Indeed, if the two first vectors are linearly 
dependent, then the third vector will always be linearly dependent. 



From Case b, it can be seen that the pair controUability/constructibihty wih 
always be guaranteed for any uniform sampUng interval. Nevertheless in the 
nonuniform case, this statement may not be true. 

More intervals for the choice of the sampling instants in systems of higher 
order can be found in jSj and yi^. According to the previous example and 
references, the underlying idea of this note is that, in spite of the restrictions 
of (17), there are intervals large enough on the time axis where the sampling 
instants can be chosen arbitrarily. Thus, this nearly free choice of the sampling 
instants can be conveniently used in conjunction with other additional criterion 
to optimize or improve different system aspects. 

4 Conclusions 

A joint characterization of reachability and observability for linear SISO nonuni- 
formly sampled discrete systems has been developed. The classical characteriza- 
tion for uniformly sampled systems appears as a simple particularization of the 
general criterion. The nonuniform sampling offers a range of situations wider 
than the uniform one in the study of these properties. 

The note stresses the nonuniform sampling sequence, which is believed to be 
an additional element for the analysis and handling of discrete systems. 
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